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ABSTRACT 

We use broadband photometry and low-resolution spectra of a complete 
sample of late-K and M dwarfs brighter than 1=22 in three helds at high 
galactic latitude to study issues relating to galactic structure and large scale 
abundance gradients in the Galaxy. The observed starcounts in each held are 
a good match to the predictions of models based on deep starcount data in 
other intermediate-latitude helds, and these models identify the late-type stars 
as members of the Galactic disk. Abundances for these late type stars are 
estimated via narrowband indices that measure the strength of the TiO and 
GaH bands in their spectra. Our results show that the average abundance in the 
Galactic disk remains close to solar even at heights of more than 2 kpc above 
the Plane. 


Subject headings: stars: late type. Galaxy: Abundances, Galaxy: Structure 


^Based partly on observations obtained at the W.M. Keck Observatory, which is operated jointly by 
the California Institute of Technology and the University of California, partly on data from the 200-inch 
telescope at Palomar mountain, owned and operated by the California Institute of Technology, and partly 
on data from the 60-inch telescope at Palomar Mountain which is jointly owned by the California Institute 
of Technology and the Carnegie Institution of Washington 
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1. Introduction 

It is by now generally agreed that the vertical distribution of stars in the Galactic 
disk cannot be represented by a simple, single-component density law. Deep, wide-held 
optical number-magnitude starcounts show that the slope of the overall density distribution 
hattens by a factor of ~ 4 between 1.2 and 1.5 kiloparsecs above the Plane, well before the 
halo becomes the dominant population at 2 : ~ 5 kpc (Gilmore & Reid, 1983; Majewski, 
1992; Reid & Majewski, 1993 - hereinafter, RM93; Gould, Flynn & Bahcall, 1996). What is 
not clear, however, is how one interprets these ‘extra’ stars, both in terms of deconvolving 
the observed distribution into thin and thicker components (Majewski, 1993; Robin, 1994) 
and in deciding the physical signihcance of the deconvolved components. 

Majewski (1993) outlines eight different scenarios for the formation of what we will refer 
to as the intermediate population (or IP II), and one can characterise the extremes as being 
those models which treat the IP II as the result of a separate, distinct star formation event 
(e.g. Wyse & Gilmore, 1988; Robin et al, 1996) and those which identify the IP II as the 
high-velocity tail of the old disk (Norris, 1987). One method of constraining these various 
hypotheses is through measurement of the abundance distribution of the stars in the IP II, 
in particular determining whether there is a signihcant gradient in the mean abundance with 
height above the Galactic Plane. To date, however, the requisite observational data are still 
relatively sparse. Yoss, Neese & Hartkopf (1987) use observations of high-latitude K-giants 
to derive a metallicity gradient of ~ —0.4 dex kpc“^ to z=700 pc. and ~ —0.18 dex kpc“^ 
thereafter (to 8 kpc), although their sample includes relatively few stars above 2.2 kpc and 
shows substantial dispersion. Trefzger, Pel and Gabi (1995) and Gilmore, Wyse & Jones 
(1995) have surveyed F and G dwarfs towards the South Galactic Pole, using VBLUW 
photometry and medium resolution, low signal-to-noise spectra respectively. Again, neither 
study include many stars above 2.5 kiloparsecs. Trefzger et al hud a relatively steep gradient 
close to the Plane (~ —0.55 dex kpc“^), and both match the K giant survey in hnding little 
variation in the mean abundance above ~ 1 kiloparsec, although Gilmore et al’s data show 
signihcantly more dispersion. Finally, Majewski’s (1992) UBV data extend to z > 10 kpc 
and suggest a somewhat higher mean abundance below 6 kiloparsecs. 

In this paper we present the first attempt to use observations of distant late-type K 
and M dwarfs to study this question. At least 80 % of the stars in the Galactic disk are 
M dwarfs, so these stars have a high surface density at faint magnitudes. Moreover M 
dwarfs can be recognized even with low dispersion spectra. At the faint magnitude level 
under consideration, any M star seen must be an M dwarf; high luminosity M giants and 
supergiants would be so far out in the galactic halo that even if they were to be found 
there in the same ratio as K giants and K supergiants (which is unlikely given the higher 
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metallicity required to produce M giants), their surface density would be insignificant 
compared to the number of bona fide M dwarfs. 

Furthermore low-abundance M-type subdwarfs can be identihed through the enhanced 
strength of the hydride bands (particularly CaH) relative to the TiO absorption (Bessell, 
1982). Until recently, the complex, diatomic-molecule dominated spectra of these stars 
dehed qnantitative abnndance analysis, bnt Gizis (1997) has shown that the atmosphere 
models developed by Allard & Hauschildt (1995) can be combined with narrowband 
spectrophotometry to derive an approximate (~ ±0.3 dex) abundance scale for these stars. 
Given these results, we have undertaken a preliminary stndy, analysing Keck LRIS data of 
a magnitude-limited sample of faint M-dwarfs. Section 2 describes the observational data; 
section 3 compares the observed colour-magnitude diagrams against Galactic starcount 
models and discusses the abundance distribntion; section 4 snmmarises our conclusions. 


2. Observations 

Our sample consists of the K and M dwarfs drawn from three fields: the 0-honr deep 
field (1=125°, b=-50°); SSA 13 (1=109°, b=74°); and SSA 22 (1=63°, b=-44°). The first 
field is 7.3 x 2 arcminntes in size (0.004 square degrees), while the other two helds are 
6' X 2' (0.003 sq. deg.). All three helds are part of a deep, K-band limited redshift survey 
being nndertaken by groups at Galtech (Gohen et al, 1996) and Hawaii (Gowie et al, 1996). 
The primary goal is a deep redshift snrvey of extragalactic objects, but objects within the 
magnitnde limits were observed irrespective of morphology. Here we concern ourselves with 
the galactic stars found during the course of these redshift surveys. This paper discusses 
observations of the complete sample of 68 stars with I <22 mag. in these helds. 

The spectroscopic sample consists of those objects with observed spectra judged to 
have sufficient signal-to-noise ratios to be useful in deducing abundances for these K and 
M dwarfs. It includes 43 stars. The spectroscopic observations of the 0 honr held sources 
are described in detail by Gohen et al (1997), while Gowie et al (1996) discuss the data for 
objects in the two Hawaii helds. 

A detailed description of the observations and analysis of the U, B, V, R, /, and K 
photometry for the 14 stars in the 0 hour held can be fonnd in Pahre et al (1997). The same 
can be found for the Hawaii helds (observed only at H, /, and K) in Gowie et al (1996). We 
have snpplemented the SSA22 data with BVI photometry obtained with LRIS on the Keck 
telescope in September, 1996, calibrating the observations using the same methods ontlined 
by Pahre et al (1997). These observations included shorter-exposure images, and hence 
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avoided saturating the brighter stars. We have used the IRAS 100 /rm maps to estimate the 
interstellar absorption in each field. Using the calibration given by Laureijs, Helou & Clark 
(1994), we estimate a total extinction of Ay ~ 0.22 mag in the 0-hour field, corresponding 
to A/ ~ 0.10, 'Ejb-i ~ 0.19 and E /_ i ^ ~ 0.08 mag; Ay ~ 0.23 mag in SSA22; and negligible 
extinction (Ay ~ 0.005 mag) in SSA13. We assume that all of the extinction lies in the 
immediate vicinity (within 200 pc.) of the Sun. Figure 1 shows the reddening-corrected 
two-colour (R — /)/(/ — K) diagram for all of the spectroscopically-confirmed stars in these 
fields. 

All of the objects were observed spectroscopically using multislit masks with LRIS (Oke 
et al, 1994) on the Keck 10-metre telescope. The slits used were 1.4 arcseconds wide , and a 
300 1/mm grating gives spectra of 16 A(~ 5 pixels) resolution covering the wavelength range 
~ 5000 to 10000 A. (The actual coverage on each spectrum depends on the exact position 
of the slit within the rectangular field.) Full details of the reduction procedures are given 
in Cowie et al (1996) and in Cohen et al (1997). The spectra from the 0 hour field were 
set on a flux scale using observations of the white dwarf standard G191-B2B (Oke 1990). 
This star was observed using a 1.5 arc-sec aperture longslit in LRIS. While the calibration 
takes out the overall spectral response of the CCD, varying losses of light at the slit lead 
to the calibrations for individual multi-slits being uncertain by up to 10 % over broad 
wavelength intervals of ~ 2000A or more. However, the widest of our narrowband indices 
spans only 250 A, so the uncertainties in the indices are dominated by the signal-to-noise 
in the individual spectra rather than by systematic errors in the flux calibration. 


3. Discussion 
3.1. Starcounts 

We have compared the observed colour-magnitude distribution of the photometric 
sample of sixty-eight stars with I<22 mag. against the predictions of starcount models. As 
the starting point for our analysis, we take model B from Reid et al (1996 - hereinafter 
RYMTS). This model includes three stellar components: the old disk, modelled as a sech^^ 
distribution, with ho of 350 parsecs; an intermediate population, with a local normalisation 
of 5 % relative to the old disk and a sech^ density law, ho = 1500 pc.; and an r“^-^ halo 
with the local density normalised to 0.15 % that of the old disk. We have also adopted a 
scalelength of 2.5 kiloparsecs for the radial density law in the Galactic disk (Robin, Creze & 
Mohan, 1992). The colour-magnitude relations adopted for these three populations are from 
photometry of nearby stars, of members of the open cluster NGC 2420 and of stars in the 
globular cluster NGC 6752, respectively. Full details are given in RYMTS, who show that 
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this model reproduces both the colour-magnitude distribution in the intermediate-latitude 
PSR 1640 held (1=41°, b=38°) and the density distribution towards the South Galactic 
Pole. 

Figure 2a compares the overall starcount predictions of model B against the 
observations in the three helds considered in this paper. The agreement is well within the 
formal Poissonian uncertainties. Limiting the sample to the magnitude range 17 < / < 22, 
we observe 19, 11 and 37 stars in the 0 hour held, SSA13 and SSA22. The predicted 
numbercounts are 15, 11 and 32 stars in the respective helds. 

Figure 2b compares the colour-magnitude distributions predicted by this model against 
the observed (I, (I-K)) distribution in each held. For greater clarity, the predicted starcounts 
in each held are matched to a solid angle of 0.01 square degrees - a factor of 2.5 greater 
than the 0-hour held and ~ 3 greater than the SSA 13 and SSA 22. The model predicts 
two distinct sequences: late F- and G-type halo subdwarfs ( (I-K) < 1.5) lying at a distance 
of ~ 10 kiloparsecs; and a sequence of late-K and M-type disk dwarfs, mainly contributed 
by the IP II, at distances of from ~ 1 to 4 kiloparsecs. The small scaleheight of the old 
disk component leads to a low surface density and an insigihcant contribution to these faint 
starcounts, even in the intermediate-latitude held SSA22. 

As described in RYMTS, the convolution of the density law and the sampling volume 
lead to a preferred distance modulus for a given stellar populaton. Since the old disk 
has the steepest vertical density law and the halo the shallowest, the former stars are 
sampled primarily at (m-M)~ 7 magnitudes, while the latter lie preferentially at (m-M)~ 
15 magnitudes. Thus, at a given apparent magnitude, the halo contributes more luminous 
(bluer) stars than the disk. The models plotted in hgure 2b clearly show that the later-type 
stars, which contribute the bulk of the stars in each held, are members of the Galactic disk. 


3.2. The abundance distribution 

The primary aim of this paper is to use the late-type dwarfs in our sample to probe 
the abundance distribution at large distances above the Plane. In order to do this we need 
to estimate both abundances and distances. We have used B-, I- and K-band data for 
individual stars within 8 parsecs to dehne (M/, (B-I)) and (M/, (I-K)) relations (hgure 3). 
The curves htted are sixth order polynomials, and the dispersion is aj = 0.35 magnitudes 
in the (B-I) calibration and aj = 0.52 magnitudes for (TK). Gomparing the distance 
estimates for individual stars, the mean diherence between the (B-I) and (I-K) photometric 
parallaxes is 40 % of the distance based on averaging the two estimates, with no evidence 
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of a systematic offset between the two distance scales. We have therefore taken the straight 
average of the two estimates as the appropriate distance for each star and adopt an 
uncertainty of ~ 30%. This is adequate for the present purposes. 

Gizis (1997) has shown that abundance estimates for late-K and M-type dwarfs can 
be derived from the relative strengths of the calcium hydride band at ~ 6800A and the 
titanium oxide 7 ( 0 , 0) band at ~ 7050A. Reid, Hawley & Gizis (1995 - RHG) have defined 
several narrowband indices designed to measure the strengths of these bands, and hgure 4 
plots data for two of the GaH indices against Ti05, the full depth of the TiO bandhead. 
In addition to data (from RHG) for single disk stars within 8 parsecs of the Sun, we have 
plotted Gizis’ measurements of the band-strength in intermediate and extreme subdwarfs. 
A comparison with the Allard and Hauschildt models shows that the former stars, plotted 
as solid square in hgure 4, have an abundance of [Fe/H] ~ —1 while the latter are 
signihcantly more metal-poor, with [Fe/H] ~ —2. We have taken nearby K and M dwarfs 
to be representative of a solar-abundance population. 

We have htted second-order polynomials to each of the three sets of stars in hgure 4 
and used these as a grid to match against narrowband indices derived from our observations 
of the faint programme stars. Table 1 lists the band-strength data for the later-type stars 
in the current sample. Broadband photometry and positions for stars in SSA 13 and SSA 
22 are given by Gowie et al (1996), while Pahre et al (1997) present UBVRIK for stars in 
the 0-hour held. Figure 5 plots the resulting bandstrengths and it is clear that the IP H 
stars in our sample lie closest to the polynomial ht to the local (solar abundance) stars, 
consistent with previous observations of IP H F and G dwarfs. We have attempted to set 
our abundance estimates onto a more quantitative scale by using the [Fe/H]=0 (nearby 
stars) and [Fe/H]=—1 (intermediate subdwarfs) sequences as hducial references. Dehning 
GaH° and GaH“^ as the values predicted by the solar and [Fe/H]=—1 relations for a given 
bandstrength, Ti05, we calculate the quantity 

^ - CaH^ 

^ ~ CaH^ - CaH-^ 

where GaH"^® is the observed strength of the GaH feature for the star in question. This 
scaling compensates for the changing separations of the two hducial sequences, although 
without an adequate number of calibrators with high-precision abundance measurements, 
we can only assume that the variation in Ac with [Fe/H] is similar throughout the full 
temperature range covered here. 

Using the distances derived from the average of the photometric parallaxes, we have 
computed the distance, z, above the Galactic Plane for each star. Limiting the sample to 
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stars with significant TiO absorption (Ti05 < 0.85), figure 6 plots the scaled CaH indices 
against z. The average values of the two indices if we include all stars are Ac 2 = —0.16 and 
Ac 3 = —0.06 (24 stars), while we derive Ac 2 = —0.27 and A^s = —0.06 for the 15 stars 
with < 2000 parsecs. As yet, we do not have sufficient calibrating observations of stars 
of known abundance to be able to transform these measured offsets into real abundances, 
but the present results are clearly suggestive of values closer to those of the old disk than 
previously suggested. 


4. Conclusions 

We have analysed photometry and low-resolution spectroscopy of faint late-K and early 
M-dwarfs in three intermediate- and high-latitude helds. A comparison with the predictions 
of a starcount model which matches the faint (R< 24) stellar colour-magnitude distribution 
in two other helds (Reid et al, 1996) shows that these stars, lying at distances of up to 2 
kiloparsecs above the Plane and members of the intermediate population. Our analysis, 
based on the strength of CaH and TiO molecular bands, demonstrates that these stars have 
abundances close to that of the old disk. Clearly, these results are preliminary, both given 
the total number of stars in our sample and the approximate nature of the calibration of 
the measured bandstrength against real abundances. However, our analysis does indicate 
a possible new direction that can be taken in exploring the abundance distribution within 
Galactic stellar populations. 


The Keck telescope project was made possible by a generous grant from the W.M. 
Keck Foundation. INR’s contribution to this research was supported partially by NSF 
grant AST-9412463; MAP was supported by NSF award AST-9157412 and the Dressier 
Foundation; DWH was supported by NSF award AST-9529170. 



Table 1. Spectroscopic indices 


Star 

Ti05 

CaH2 

CaH3 

PDOK 03 

0 hour field 

0.765 

0.604 

0.816 

PDOK 04 

0.792 

0.627 

0.837 

PDOK 09 

0.550 

0.494 

0.763 

PDOK 16 

0.643 

0.503 

0.790 

PDOK 19 

0.597 

0.597 

0.811 

PDOK 21 

0.458 

0.428 

0.726 

PDOK 22 

0.549 

0.444 

0.712 

PDOK 49 

0.387 

0.419 

0.693 

PDOK 50 

0.625 

0.626 

1.003 

PDOK 91 

0.549 

0.479 

0.718 

PDOK 97 

0.333 

0.348 

0.594 

PDOK 101 

0.392 

0.597 

0.689 

17 

SSA 13 

0.461 

0.433 

0.697 

35 

0.435 

0.347 

0.655 

73 

0.417 

0.508 

0.808 

84 

0.378 

0.321 

0.583 

97 

0.583 

0.466 

0.709 

24 

SSA 22 

0.290 

0.333 

0.387 

61 

0.644 

0.400 

0.845 

112 

0.400 

0.475 

0.754 

157 

0.489 

0.477 

0.755 

169 

0.615 

0.435 

0.713 

171 

0.601 

0.436 

0.706 

200 

0.616 

0.480 

0.779 
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Fig. 1.— The (B-I)/(I-K) two-colour diagram for stars in the 0-hour field (solid dots), SSA 
13 (crosses) and SSA 22 (solid triangles). The reference sequence is provided by Leggett’s 
(1992) photometry of stars within 8-parsecs of the Sun (plotted as open squares). The 
error-bar shows representative uncertainties for I~ 21 magnitude and fainter. 
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Fig. 2.— (a) Predicted and observed number counts for the three fields in this study. The 
solid histogram outlines the total counts; the dotted line marks the contribution of the old 
disk; and the dashed line separates the contribution from the IP II and halo. 
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Fig. 2.— (b) Predicted and observed (I, I-K) colour magnitude diagrams. The right-hand 
panel show the predictions for each field (and for a solid angle of 0.01 square degrees) made 
by the model described in the text, where the red sequence is due to old disk (crosses) and 
IP II stars (triangles) and the blue due to halo stars(open circles). The left-hand panels plot 
the observed distributions. 













Fig. 4.— Calcium hydride strength as a function of the TiO 7 ( 0 , 0) bandstrength. Data for 
single stars within 8 parsecs of the Sun (from RHG) are plotted as open squares. These dehne 
a solar-abundance sequence. The solid symbols are subdwarfs from Gizis (1996), with the 
squares being ‘normal’ subdwarfs with [Fe/H] ~ —1 and the triangles ‘extreme’ subdwarfs, 
[Fe/Hl r\j — 2. Second-order polynomials htted to each sequence are also shown. 















